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Abstract 

We argue that the paradoxal softness of the red blood cells (RBC) in fluctua- 
tion spectra experiments is apparent. We show that the effective surface shear 
modulus jx s of the RBC obtained from fluctuation data and that measured 
in static deformation experiments have the same order of magnitude. A sim- 
ple micromechanical model of the RBC developped for this purpose accounts 
for the influence of a finite-thickness cytoskeleton on the fluctuations of the 
composite membrane-cytoskeleton system. The spectrin network cytoskele- 
ton with the bulk shear modulus estimated as \x ~ 105 -j- 165 Pa contributes 
to both normal and tangent fluctuations of the system and confines the fluc- 
tuations of the lipid membrane. The ratio of mean square amplitudes of the 
RBC normal and tangent fluctuations < X 2 n > / < X% > calculated in the 
frame of the model is 2-3 orders of magnitude smaller that it is in the free 

membrane with the same bending and shear moduli. 
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The mechanical properties and fluctuation spectrum of the living cell depend strongly 
on the coupling between its fluid phospholipid membrane (ME) and stiff cytoskeleton (CS) 
composed of the cross-linked biopolimer networks [1-4]. In red blood cells (RBC), investi- 
gated intensively by physicists since many years [5-7], the coupling between the lipid bilayer 
and the spectrin filaments network provides the ME with viscoelastic properties. In addi- 
tion to its bending modulus K ss 1CT 19 — 1CT 20 J, the RBC ME acquires two dimensional 
(2D) static shear modulus revealed by micropipette suction [5] or by other types of static 
deformation response. Measurements of thermally activated shape fluctuations [8] refine 
the picture and give the access to the dynamic properties of RBC. The basic theoretical 
model [9] explaining both statics and dynamics of RBC considered the ME as practically 
independent on the CS influence, but with strongly variable 2D shear modulus. Namely, 
thermally excited surface undulations of RBC (the cell flickering) [8] was described as a 
motion of a fluid ME with /j, s rs 0. Thus the free ME model [9] has established the opinion 
that RBC combine paradoxally liquid-like and solid-like mechanical properties, with solid- 
like behavior being manifested only during large-scale shape changes [9-10]. However, the 
modern microrheological technique [e.g. 11-12] has shown that mechanical behavior of RBC 
satisfies Hook's law in a broad region of deformation [13]. This result is incompatible with 
both deformation scale separation principle and non-linear shear modulus behavior resulting 
from the model [9]. 

As it will be shown later apparent softness of RBC obtained from the fit of thermal fluc- 
tuation spectra by the free ME model is related to the simplified form of its elastic energy. 
In [9] it is reduced to the ME bending and shear terms only. In recent theoretical works 
[14-17] it was proposed to take into account the influence of viscoelastic CS properties on 
the coarse-grained mechanics of the composite ME-CS system. Energy of surface tension 
and especially the term responsible for the ME confinement were introduced in the system's 
Hamiltonian [16-17]. Resulting model clarified the effective bending modulus behavior dur- 
ing normal fluctuations. Besides, it demonstrated the physical factors which limit normal 
fluctuations of the ME. By contrast, model [16-17] could not explain the main RBC paradox 
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: the apparent change of its shear properties from the solid-like to the liquid-like ones. In 
approach [16-17] the CS was treated as an infinitely rigid shell (i.e. with an infinite shear 
modulus) making impossible the calculation of the effective shear modulus for the ME-CS 
structure. Infinite rigidity of the CS prevents any tangent fluctuation in the model, and 
thus does not allow to compare its predictions with the results of the free ME model [9]. 
Note that the conclusion about liquid-like behavior of RBC was obtained in [9] from the 
calculation of the ratio < X\ > / < Xj > of mean-square amplitudes of normal and tangent 
fluctuations (see Fig. 1). To solve the problem authors of Ref. [18] proposed to consider the 
dependence of /i s value on the wave vector. In their model /i s modulus tends to zero only 
for the fluctuations with wavelengths smaller than the mesh size of the spectrin network. 
Nevertheless, the < X\ > / < Xj > ratio has not been calculated. 

In the present work we propose more realistic continuum model of the composite ME- 
CS system which brings a solution to the paradoxal problem of the RBC apparent softness 
in thermal fluctuation measurements. We show that the 2D shear modulus of the RBC 
obtained from the fluctuation data and that measured in static deformation experiments 
have the same order of magnitude /x s 6 x 1CT 6 N/m. For that purpose, in addition to the 
ME confinement, we take into account the contribution of normal and tangent fluctuation 
modes of a finite-thickness CS to the surface free energy of the ME-CS system. These 
modes, determined in the simplest case of incompressible CS by its bulk shear modulus /x 
only, are rapidly attenuated in the cell interior. Normal CS fluctuations modify the ME 
bending properties and tangent fluctuations supply the composite system with a finite 2D 
shear modulus /i s . Mean-square amplitudes of normal and tangent fluctuations of the system 
are calculated in the framework of the proposed model. They are perfectly consistent with 
the RBC flicker spectroscopy data [8,19] but their ratio < X\ > / < Xf > is 2-3 orders 
of magnitude smaller than that attributed in [9] to the solid-like ME. This result shows 
that the RBC fluctuational softness is apparent and its shear properties in both static and 
dynamic experiments are the same. 

The material constants for these estimations were obtained by a double-step numerical 
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fit of the effective bending rigidity dependence on the wave vector q [19]. The same type of 
experimental data is frequently used in the theoretical estimations of the RBC mechanical 
properties [16,17,20]. The data being obtained using classical Fourier transform, the models 
are usually formulated in the framework of plane geometry. We also follow this approach but 
limit the region of its applicability. A model which does not take into account the curvature 
of the system is valid [14-17] only for a treatment of fluctuations with the wavelengths 
X « R, where R « 4 fim is the effective RBC radius. The description of the ME-CS 
system motion is performed then in terms of continuum mechanics. This approximation is 
favored by the fact that the attachment points between the spectrin network and the bilayer 
are rather rare and their area is negligible (Fig. 1). Besides the actin nodes there exist 
additional ME-CS connections through ankryn complexes distributed randomly. Together 
with defects and dynamic structural rearrangements in the spectrin network they contribute 
to the entropic elasticity of the RBC CS and justify the continuum description. Finally, 
we consider the fluctuation wavelength region for which the influence of active processes is 
reduced to the increase in the effective temperature [21-22]. 

The free energy of normal thermal fluctuations of a ME-CS system contains the fluctu- 
ation energy of a flat ME and the CS contribution W: 



Here u m is a normal displacement of the ME, K, a and 7 are the coefficients of bending, 
surface tension and pinning of the ME, respectively; u c is a displacement field of the CS. 

The ME fluctuation energy (first three terms in Eq. (1)) is limited to the sum of bending, 
surface tension and ME confinement energy contributions. In a plane geometry adopted in 
the present work bending fluctuations are independent on stretching and the energy due to 
the bending-stretching coupling can be omitted, though in a general system with curvature 
it cannot be neglected [15,23,24]. The CS is considered as an incompressible viscoelastic 
plate (Fig. 1) with the thickness h ~ 0.03 fim and with a finite 3D shear modulus //. The 
last term in (1) expresses the surface density of the CS fluctuation energy: 




(1) 
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W = jf dzfie^, (2) 

where = ^(l^ 1 + l^ 2 -) is the CS strain tensor, u\ = u x ,u y ,u c z and = x,y,z; and z is 
directed perpendicular to the ME-CS interface. Due to the small Reynolds number of the 
system the motion of both ME and CS is overdamped. Thus the state of CS is determined 
from the motion equation of an incompressible medium with neglected inertia terms: 

Vp = /iAu c + r]Au c (3) 

with p standing for pressure, rj for the CS viscosity and u c for the CS velocity field. In 
a more general case of a compressible isotropic medium the energy (2) would have two 
terms related to two independent elastic Lame coefficients A and /i. The equation of motion 
of a compressible isotropic medium would contain the term A grad div u, where u is a 
displacement field. This term stands for the density of forces of isotropic deformation. In 
the incompressible medium div u = and A tends to infinity. As it is shown in [25] the limit 
value of the term A grad div u in the incompressible medium is equal to -Vp. This leads to 
the motion equation of the incompressible medium in the form given by Eq. (3). 

Thermally excited CS fluctuation modes are rapidly attenuated in the direction perpen- 
dicular to the ME plane. The solution of Eq. (2) compatible with bending fluctuations of 
the ME and satisfying div u c = condition describes the mode of this type: 

u c x = —iU q qz exp(—qz + iqx); u c z = U q (l + qz)exp(—qz + iqx); u c y = 0; (4) 
p = —2iqU q (r]u + i/i)exp(—qz + iqx); (5) 
u c = ~icou c . (6) 

The wave vector of fluctuations in (4-6) is chosen along the x axis, U q expresses the am- 
plitude of normal displacements at the ME-CS interface z = (see Fig. 2). All tangent 
displacements in the mode of this type vanish at the interface. 

In the present work we consider the physical situation, where at the interface the CS 
velocities u c and displacements u c coincide with those of the ME: u c | z=0 = du m /dt and 
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u c | 2=0 = u m . In a more complicated case of ME-CS mutual gliding, only the normal com- 
ponents of both velocity and displacement fields are preserved at the interface and a viscous 
drag influences tangent fluctuation modes. However, in the plane geometry approximation 
this does not affect the normal fluctuation modes (which are of primary interest in the RBC 
model) and maintains the relation u m = U, where u m and U are the normal displacement 
of ME and CS, respectively. 

Finally, the contribution W q of the CS to surface fluctuation energy (1) is obtained by 
substitution of bending displacement field (4) into elastic energy (2) and following integration 
over the CS thickness: 

W q = fiq\U q \ 2 (l - exp(-2qh)(2q 2 h 2 + 2qh + 1)) . (7) 

After development of free energy functional (1) into Fourier series we get the total energy 
of normal ME-CS fluctuations in the form: 

E = SY,E(q)\U q \ 2 (8) 
q 

suitable for analysis of the RBC fluctuation spectrum [8,19]. Function E(q) is expressed as 

E(q) = Kq 4 + aq 2 + 2W q /\U q \ 2 + 1 . (9) 

The sum in Eq. (8) runs over all possible q in a square plate with the side equal to nR; 
the interface area S is taken to be (irR) 2 . Note that in the limit qh — > CS fluctuations 
simply renormalize the ME bending rigidity: W q ~ 4/x/i 3 |[/ g | 2 g 4 /3; and in the limit qh — > oo 
CS fluctuation wave transforms into the wave of normal fluctuations on the surface of a 
semi-infinite medium: W q ~ jj,q\U q \ 2 . Due to a strong exponential attenuation of solution 
(4-6) with the depth z the last approximation is suitable in a wider region q > h. 

The authors of [19] give the fluctuation spectra of RBC as the dependence of — log < 
\U q \ 2 > q A on the wave vector q in the wavelength region from 0.25 \xm to 2 jim. They stress 
also that the region of fidelity of their results does not exceed 0.5 fiin < A < 1 /jm limits. 
Long-wavelength limit is determined by the validity of the plane geometry approximation 
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(used in [19] for the data processing). Lack of fidelity in the short-wavelength region is 
related to the experimental technique limitations. Therefore, we performed the fit of the 
RBC fluctuation spectrum in the fidelity region only. Using equipartion E(q) oc< \U q \ 2 > _1 
one can fit the expression q 4 /E(q) proportional to the experimental value q 4 < \U q \ 2 >, 
and usually related to the effective bending rigidity [16,17,20]. Double-step least-square 
fit (Fig. 3) of relative values of q 4 < \U q \ 2 > leads first to the expressions of the RBC 
material constants in function of bending modulus K. The inset in Fig. 3 shows the relative 
contributions of the bending, surface tension, CS fluctuation and confinement energies to 
the total mode energy E(q) for one of the three cells studied in [19]. Similar relation 
between different contributions exists for two other cells. Then, taking into account the 
average value K = 3.4 x 10~ 20 J estimated in [19] we obtain the absolute values of the 
constants: « 4.9 x 1(T 20 J; crW « -5.9 x 1(T 6 N/m; 7 « « 3.1 x 10 8 N/m 3 ; ^ « 148 
N/m 2 ; « 2.2 x 10~ 20 J; « -2.3 x 1(T 6 N/m; 7^ « 1.8 x 10 s N/m 3 ; ^ « 105 
N/m 2 . The fit of experimental data [19] for the third cell (not shown in Fig. 3) results in 
/i® 165 N/m 2 . The found value of the CS 3D shear modulus /i is in a good agreement 
with an effective spring constant of spectrin filament ~ 4 x 10~ 6 N/m in 30 nm thick CS. 
The effective value of the 2D shear modulus /i s can be estimated in the simplest way as 
fi a ^fjih= (3.2 + 5) x 10- 6 N/m. 

Let us now show that to explain the RBC behavior in thermal fluctuation experiments 
we do not need to suppose vanishing of its effective 2D shear modulus (// s « 0). To do 
this we reexamine two main arguments of the free ME model [9] which have lead to the 
conclusion about negligible /i s value and, then reanalyze the same experimental data in the 
framework of the composite ME-CS model developed in the present work. 

The arguments developed in [9] can be resumed qualitatively in the following way: 1) 
The ratio of mean-square amplitudes of normal and tangent fluctuations < X 2 > / < X 2 > 
in the model [9] has the order of magnitude ~ [i s R 2 /K. 2) Thermal fluctuation experiments 
show that the RBC thickness fluctuation profile has the peak near the rim of the cell [8] . Due 
to the RBC shape the main contribution to the thickness fluctuations near the rim comes 
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from tangent fluctuations (see Fig. f). Thus, in a system with ji s pa and, consequently, 
with strong tangent fluctuations, the peak should find itself near the rim. With the value of 
ji s pa 6 x 10~ 6 N/m known from static experiments the model [9] gives the < X\ > / < X? > 
ratio 2-3 orders of magnitude greater than the experimental one. The discrepancy is then 
attributed to the high value of /x s and to preserve the model the effective 2D shear modulus 
is considered to be vanishing. 

Estimation of the < X\ > / < Xf > ratio in the framework of the present model 
gives however quite different result which is consistent with the static /x s value. The main 
difference with respect to the model [9] is determined by the contribution to total energy 
(8) from the CS, both from its fluctuations and from the ME confinement. Mean-square 
amplitudes of normal and tangent fluctuations of the composite ME-CS system can be, in 
principle, calculated using equipartition theorem. To avoid however, a cumbersome discus- 
sion of different mode polarizations for tangent fluctuations we prefer to illustrate the results 
in the way proposed in [26] and based on the fluctuation-dissipation theorem. Mean-square 
amplitudes of normal and tangent fluctuations are then determined by static normal and 
tangent response functions, respectively [26,27]: 



where response aii(u) defines amplitude Ai(u) of ME particle motion under periodic external 
force Fi = F®exp(—iut) application: = Ai/Ff. The values of responses at zero frequency, 
and consequently, the amplitude of mean-square fluctuations are independent on the system 
dissipative properties. The responses of a flat ME have been discussed in detail in [14,15]. 
Using the same formalism we obtain normal response a n (0) of the composite ME-CS system 
in the form: 



where E(q) is defined by Eq. (9). The minimal wave vector q min of fluctuation which 
contributes to response (11) is determined by the RBC finite size: q m in = 1/-R- Upper cut- 
off q max is defined [27] by the fact that the number of fluctuation modes in the system is finite. 



< X 4 (0)X,(0) >= K B Tai(0), 



(10) 




(11) 
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In the lipid ME it is determined by the total number N of lipid molecules: q max = vNq min . 
In such a case the number of normal fluctuation modes is equal to N. 

The expression for the tangent response a t (0) is obtained along the same line. In the 
general case the function ct t (u) can be presented as a t (uo) = ct s t h {u) + ctf(u), where a s t h (uo) 
and af{uj) are contributions of the shear and stretching fluctuation modes, respectively [14]. 
However, the lipid ME is commonly considered as incompressible. In such a case of{ui) = 
and a t {uj) = a s t h {u). Note, that shear fluctuations in the proposed model depend mainly 
on the CS properties. To obtain the energy of shear modes we use the solution of Eq. (3) 
compatible with shear fluctuations of the ME: 



u 



x 



0; Uy = U q exp(—qz + iqx); u c z = 0; p = 0; u c = — iuju c . (12) 



Similar to mode (4-6), thermally excited shear CS fluctuations are rapidly attenuated in the 
^-direction (see Fig. 2). Substitution of displacement field (12) into elastic energy (2) and 
integration over the CS thickness results in the shear fluctuation energy 

E sh = Sj2w\U q \ 2 (l-exp(-2qh)) . (13) 
q 

Corresponding shear contribution to the tangent response of the composite system is then 
expressed as: 

a sh (0) = / — (14) 

* J qmin 47r//(l - exp(-2qh)) ' 

To calculate responses (11) and (14) we use the RBC material constants K, a, 7 and fi 
fitted above and take iV xs AA x 10 7 molecules [28]. The normal response a n (0) is estimated 
to be ~ (6 4- 9) x 10 4 m/N. This value is 2-3 orders of magnitude smaller than the response 
function of the free ME model with the same bending rigidity. For the tangent response 
function a t (0) we obtain a t (0) ~ (8 4 11) x 10 5 m/N. Thus, final < X 2 n > / < X 2 t > 
ratio is also 2-3 orders of magnitude smaller than that in the free ME model and has no 
contradiction with the thickness fluctuation profile experiment. Let us stress, that the value 
obtained in the approximation of an incompressible ME-CS system is an upper limit for 



the < X^ > I < X? > ratio. The system compressibility increases the < X t 2 > value and 
preserves the < X\ > one. It means that for a compressible system our approach works 
even better. 

Knowledge of responses (11) and (14) allows us to make a little better estimation of 
fi s modulus than that given by the relation /i s jih. The following estimation is more 
appropriate for the static microrheological experiments using the probe particles technique 
(see, for example [13]). In such experiments the cell is usually replaced for data processing by 
a 2D solid-like ME. The static shear properties of the both systems are roughly equivalent 
provided a$ h (0) = P^ h (0), where ffi h (Q) denotes the shear contributions to the tangent 
response of the 2D solid-like ME. The response Pt h (co) of the 2D ME in connection with the 
PP microrheology was studied in [14,15]. Its value 



depends on the wave-vector cut-off qpp = 2n/Rpp [14,15], where Rpp is the radius of the 
contact area between the RBC and the PP (i.e. silica bead). The equality between a^ h (0) 
(see Eq. (14)) and Pf h (0) results in the following effective 2D shear modulus // s : 



Taking R PP ps R/20 we obtain that the CS with the bulk shear modulus /i ps 105 + 165 Pa 
and h « 30 nm thickness induces in the composite ME-CS system the 2D shear modulus 
[/, s ~ (5.1 -r 8.1) x 10~ 6 N/m. This value is very close to that obtained in static deformation 
experiments [5-7]. 

In addition to the resolution of the RBC apparent softness paradox, the composite ME- 
CS model has a striking feature: negative a value obtained from the fit of fluctuation spectra 
[19]. Negative osmotic pressure difference AP — 2a /R 1 1.5 Pa insures a good contact 
between the ME and the CS. It makes the mechanical interaction between the parts of the 
composite system much less dependent on the attachment points. In such a case the CS 
induces (through the controlled tension mechanism) additional ME fluctuations which can 




(15) 




(16) 
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even result in an elastic instability of the cell [29]. A high negative surface tension would 
break the symmetry from the normal biconcave RBC to the echinocyte shape. On the other 
hand, a high positive surface tension corresponding to a sufficiently big positive osmotic 
pressure difference would make a normal cell unstable with respect to the spherical shape. 
These points attract considerable attention in the cell biology field [30,31]. 

More detailed physical discussion of the surface tension could be done in the frame of 
a more rigorous model taking into account both finite shear modulus /i s and the real RBC 
shape. In such a model surface tension becomes a tensor characteristic (y^ of the ME that 
satisfies the Laplace equilibrium equation: agg/Rg + a^/R^ = AP. Here Rg and R^ are the 
RBC curvature radii along 9 and spherical coordinates, respectively. Due to the rotational 
symmetry of the non-pathological RBC shape <7g^ = 0. The values of the ogg and 
components are dependent on the 9 coordinate of the ME point. Thus, change in signs of 
principal curvatures 1/Rg and l/R^ can alternate compressed and stretched regions at the 
cell surface. If AP is negative (as it is obtained in the present model) then the ME region 
near the rim is compressed: positive curvatures near the rim lead to <jgg < and < 0. 
This fact favors normal fluctuations and increases their contribution to < 5d 2 > (see Fig. 1). 
Near the cell center the RBC shape is locally concave and both curvatures are negative. The 
ME is stretched in this region (agg > and > 0) and its normal fluctuations are reduced. 
Consequently, the amplitude of normal fluctuations increases with the distance from the cell 
center. Remarkably, it puts negative osmotic pressure value in a good agreement with the 
location of the thickness fluctuation profile peak. 

In conclusion, we showed that the RBC paradox, consisting in solid-like elastic behavior in 
static suction and liquid-like behavior in thermal fluctuations experiments is apparent. The 
influence of the finite-thickness CS which confines the ME and contributes to the fluctuation 
energy of the composite ME-CS system reduces considerably normal fluctuations of the RBC. 
Contrary to the predictions of the free ME model, the amplitude of normal fluctuations in the 
composite system is smaller than the amplitude of tangent fluctuations. Resulting effective 
2D shear modulus of the RBC shows in thermal fluctuation spectra the value typical for the 
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solid-like elastic behavior. 

The authors thank G. Mennessier for helpful discussions. S.B.R. is grateful to the Bian- 
nual Program of the French Minestere de l'Education Nationale, de l'Enseignement Superieur 
et de la Recherche for financial support. 

Figure captions 

Fig. I. Equilibrium shape and fluctuations of the RBC. Mean-square amplitude of the 
RBC thickness fluctuations < 8d 2 > depends on the mean-square amplitudes of normal 
< X 2 > and tangent < X 2 > fluctuations and on the cell surface point: < Sd 2 >=< X 2 > 
cos 2 (3+ < X 2 > sin 2 (3; here (3 is the angle between the vertical direction and the RBC 
surface normal. Near the rim of the cell the main contribution to < 5d 2 > comes from 
tangent fluctuations. Inset: Schematic representation of the RBC membrane and finite- 
thickness cytoskeleton. 

Fig. 2. Displacement fields in the cytoskeleton. 2D sections by the plane containing z 
and x axes are shown. The three different modes presented are characterized by the same 
wave vector q. The membrane is located in (x, y) plane at z = level. According to the 
boundary condition (see in the text), the displacement of the membrane surface is equal to 
that of the cytoskeleton at z = level, therefore the membrane is not shown, (a) Bending 
mode, (b) Shear mode. The displacement field is perpendicular to the plane of the figure. 
Two opposite directions of the displacement field are shown by crosses and full circles. Their 
size is proportional to the displacement value. 

Fig. 3. Wavelength dependence of q 4 /E(q) (E(q) is given by Eq. (9)) related to the 
RBC effective bending rigidity in the flicker spectroscopy experiments [19]. Fits (solid lines) 
of the data for two different cells (rhombuses and squares) are presented. Inset: Calculated 
wavelength dependence of the terms in the normal fluctuation mode energy E(q): bending 
contribution E\ = Kq A , surface tension contribution E 2 = crq 2 , CS fluctuations contribution 
E 3 = 2W q /\U q \ 2 and confinement E± = 7. Though the surface tension term is negative, the 
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globally positive mode energy E(q) insures stability of the system with respect to fluctua- 
tions. 
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